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TEST-STAND STUDY OF A VEHICLE UPPER STAGE EMPLOYING A 


ROTATING-SO LID-ROCKET ATTITUDE CONTROL SYSTEM 

AND A COMPARISON WITH ANALOG STUDIES 

By A. Thomas Young and Robert L. James, Jr. 
Langley Research Center 


SUMMARY 


A dynamic -test -stand investigation in three degrees of rotational freedom 
of a space-vehicle upper stage employing a rotating- solid-rocket attitude con- 
trol system is described. Experimental data received from an onboard telemetry 
system and a unique photographic technique are presented and compared with 
results obtained from an analog simulation of the dynamical system. Analysis 
and comparison of the results from the various sources showed that the system 
performed as predicted by the analog simulation, thus validating this mathe- 
matical representation of the dynamical system. It is shown that the telemetry 
results, although limited in accuracy, are very valuable in analyzing the atti- 
tude control characteristics of the space-vehicle upper stage. 


INTRODUCTION 


The rotating- solid-rocket control system which is the subject of investi- 
gation in the present report has certain attractive features which more con- 
ventional spacecraft control systems do not possess. Jet vanes, jetavators, 
movable main rocket nozzles, and secondary fluid injection systems must be 
developed as an integral part of the main rocket. However, because of the size 
and complexity of most propulsion systems, it is desirable to have a control 
system which can be independently developed. The present system which utilizes 
four auxiliary rockets for control forces can be developed independently of the 
main rocket, and in addition, it is inherently capable of providing velocity 
control, retrothrust, and spin-up forces. This control concept was developed 
under NASA contract for the George C. Marshall Space Flight Center (see ref. l) , 
and the Langley Research Center is conducting a comprehensive investigation of 
this promising control- system concept. 

An analog study of the application of this system to control the attitude 
of a space-vehicle upper stage is presented in reference 2. In this analog 
analysis, mathematical expressions were used to represent the physical systems 
comprising the space-vehicle upper stage. The present investigation was under- 
taken to determine the adequacy of this analog representation of the space- 
vehicle upper stage. 


I 


A full-scale model of the space-vehicle upper stage which included the 
electro -mechanical flight systems was mounted in a dynamic test stand which 
provided three degrees of rotational freedom. The attitude control system was 
activated and provided control during programed vehicle maneuvers and disturb- 
ances imparted to the vehicle by solid- fuel disturbance rockets. Data received 
from an onboard telemetry system provided information on control-rocket deflec- 
tions and guidance signals from attitude and rate gyros. In addition, the 
control-rocket thrust performance was monitored with the use of pressure trans- 
ducers. A unique camera technique was employed to obtain the Euler angles 
which defined the vehicle attitude during the dynamic test. The data obtained 
from this experiment were compared with the results obtained from an analog 
simulation of the dynamic model in order to evaluate the adequacy of the math- 
ematical analysis, and these comparisons are reported herein. 


SYMBOLS 


A,B,C triangle defined by pivot point of dynamic model, position- 

indicating ball, and camera 

e distance between graduated framework and position- indicating ball 

when dynamic model is vertical, in. 

F thrust of one control rocket multiplied by cos 10° (control 

rocket cant angle), lb 

h distance between camera and graduated framework, in. 

I X ,I Y ,I Z mass moment of inertia about X-, Y-, and Z-axes, respectively, 
slug-ft^ 


M x ,My,M z 


disturbance moments about X-, Y-, and Z-axes, respectively, ft-lb 


0 

P 

P 

r 


t 

X,Y,Z 


origin of Xe~, Ye-, and ZE-axis system 

projection of position-indicating ball into y e ,z e plane 

rolling velocity, deg/sec or rad/sec 
pitching velocity, deg/sec or rad/sec 
yawing velocity, deg/sec or rad/sec 

distance between pivot point of dynamic model and position- 
indicating ball, in. 

time, sec 

body axes of dynamic model 
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Xe,Ye.,Ze earth-fixed axes 


Xj>,yg,z-g distance along earth-fixed axes, in. 

Xr,Yr,Z;r reference axes 

x cr center of rotation from nose, ft 

Xp moment arm for control moment about X-axis, ft 

Xq moment arm for control moment about Y-axis, ft 

x r moment arm for control moment about Z-axis, ft 


a,p,o> angles contained in triangle defined by A, B, and C, deg 

8 control- rocket deflection angle, deg 

T) angle between camera line of sight and longitudinal axis of 

dynamic model, deg 

pitch, yaw, and roll Euler angles of dynamic model relative to 
earth-fixed axes, deg 

9 R> 1 J r R>0R pitch, yaw, and roll reference angles relative to earth-fixed 
axes, deg 

pitch, yaw, and roll attitude error angles relative to reference 
axes, deg 

Subscript : 


o 


initial value 


Dots over symbols denote differentiation with respect to time. Primes 
indicate uncorrected photographic data. The numbers 1, 2, 5, and A denote the 
four control rockets as shown in the following sketch (view looking forward 
along the X-axis): 
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MODEL AND TEST APPARATUS 


General 

A dynamic model consisting of a spacecraft and ballast was used in the 
test- stand investigation to simulate a typical space- vehicle upper stage. (See 
fig. 1. ) The spacecraft housed the control system of research interest, a 
guidance system, a telemetry system, and other supporting units. A detailed 
description of the guidance and control systems is given in reference 2. The 
ballast was employed to provide a dynamic model with center-of-gravity and 
inertia characteristics that approximately represented the space-vehicle upper 
stage used as a test bed in reference 2. The study was conducted with the use 
of a dynamic test stand (see fig. l) that allowed three degrees of rotational 
freedom. The test setup consisted of the dynamic model mounted on a hydraulic 
bearing pivot allowing three degrees of rotational freedom relative to and 
within the bounds defined by the fixed test-stand structure. The desired test 
data were obtained by the telemetry system and a unique camera technique. 


Dynamic Model 

A sketch of the dynamic model and fixed test-stand structure is included 
in figure 2, and a block diagram showing the integration of the individual sys- 
tems that form the dynamic model is given in figure 3- 

The point of rotational freedom was located at station 132 inches which 
was also the approximate location of the center of gravity of the dynamic model. 
This geometric arrangement was accomplished by adding an additional constraint 
on the ballast mass and position selection. The criterion for the ballast 
design was to satisfy the center-of-gravity requirement as well as to produce 
pitch and yaw inertias similar to those of the space-vehicle upper stage. 

The pitch and yaw inertias were each 1065 slug-ft^, and the roll inertia 
was 55.4 slug-ft^. These inertias were nearly constant since the mass change 
during the test was small compared with the mass of the dynamic model. 

Control system .- The control system to be considered herein is pictorially 
shown in figure 4. The thrust of four end-burning solid-propellant rocket 
motors provides the control forces. The control rockets are bearing mounted 
in a minimum-volume configuration so that rotation of two motors produces pitch 
control moments, rotation of the remaining two motors produces yaw control 
moments, and differential deflection of all four control rockets produces roll 
control moments. Figure 5 schematically shows the control rockets deflected 
to correct a positive yaw error (top view) and a positive pitch error (side 
view). The control-rocket deflections were limited by design to ±60°. The 
torque to produce a control-rocket deflection is provided by a direct-current 
electrical drive motor. The control system has four drive motors, one geared 
to each control rocket. Also geared to each control rocket is a three-deck 
potentiometer. One deck provides control-rocket position feedback for stabil- 
ity, and the remaining two are used as data sensors to provide information 
concerning the control-rocket positions. 
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Guidance system . - A position and rate monitoring proportional- guidance 
system is employed in the experimental analysis of the control system previously 
discussed. The guidance system uses two two -degree-of -freedom gyros which pro- 
vide inertial attitude information in roll, pitch, and yaw relative to a ref- 
erence axis system defined hy the spin vectors of the gyros. The guidance sys- 
tem also employs three rate gyros which are fixed in the spacecraft such that 
they measure "body rates. The outputs of each gyro and the control-rocket posi- 
tion feedbacks are weighted and appropriately summed to provide error signals 
to the control system. (See fig. 3 ») The numerical values and procedure for 
the weighting and summing are described in reference 2. 

Telemetry system .- The purpose of the telemetry system was to monitor the 
desired performance parameters and telemeter these data to a receiving sta- 
tion. The system was typical of units employed in sounding-rocket flight-test 
research. The rate and attitude information was obtained from gyros contained 
in the guidance system. 

The results which were obtained by the telemetry system are 

(1) Roll, pitch, and yaw body rates from rate gyros 

( 2 ) Roll, pitch, and yaw attitude error angles from attitude gyros 

(3) Rotational position of each control rocket 

( 4 ) Chamber pressure of control rockets 


Test Apparatus 

Test stand .- The test stand used in the experimental study is shown picto- 
rially in figure 1 and schematically in figure 2 . It provided a fixed rigid 
pivot about which the dynamic model was free to rotate. A ball and socket into 
which hydraulic fluid was forced under pressure constituted the pivot point. 
Since the hydraulic-fluid viscosity was the only source of resistance, the 
pivot-point friction was negligible. The point of rotational freedom was held 
fixed by structure securely mated with a concrete base. 

Camera . - Surrounding the test stand was a structure which provided support 
for the photographic apparatus. (See fig. 1 .) A l6-mm camera which operated 
at a rate of 64 frames per second was mounted above the dynamic model and viewed 
the model from above looking down. Also supported by the structure was a grad- 
uated framework mounted between the camera and the nose of the dynamic model. 

Disturbance rockets . - Solid- fuel rocket motors with an average thrust level 
of 3b pounds each and a thrust duration of about 10 seconds were mounted on the 
dynamic model to produce known disturbances. Two disturbance rockets were 
located near the nose at station 36.5 inches, and two were mounted on the lower 
portion of the ballast at station 163 inches. The upper rockets were alined 
90° apart, and the lower rockets were alined l80° apart with a distance of 
10 inches between the center lines. 
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Restraining explosive "bolts .- The dynamic model was restrained in a verti- 
cal position at the initiation of the experimental testing by two explosive 
bolts. These bolts were attached between the base of the dynamic model and the 
concrete pad located beneath the test stand. The dynamic model was free to 
rotate after pyrotechnic severing of the two restraining bolts. 


TEST DESCRIPTION 


The sequence of events for the experimental test is shown in the following 
table : 

Time, 

sec 


Precess attitude gyro spin axis to define reference axis system -^-A 

Zero time for the programer used in the test 0 

Control-rocket ignition 76.5 

Dynamic -model release 78.5 

Initiate roll attitude and rate control as well as pitch and yaw 

rate control 79-5 

Initiate pitch and yaw attitude control to give full three-axis 

control 93.5 

Ignite yaw disturbance rocket 99*5 

Ignite roll disturbance rockets 112.9 

Ignite pitch disturbance rocket 126. 5 


The initial orientation of the dynamic model was vertical and it was held in 
this position until approximately 2 seconds after control-rocket ignition at 
which time it was released by pyrotechnically severing the two restraining 
explosive bolts. After release the model was floating free with the only sta- 
bilizing forces being produced by the research control system. 

The dynamic model is controlled from 0.8 second after release to termina- 
tion of the experiment. The control period is divided into two phases. The 
first phase is defined as roll control and is effective for the first 14.2 sec- 
onds. During this time the model is subject to roll attitude and rate control 
as well as pitch and yaw rate control. The second phase, which is effective 
from t = 93*5 seconds to the termination of the experiment, is a three-axis 
control phase. During the second phase the model is controlled in attitude and 
rate about the roll, pitch, and yaw axes. A detailed discussion defining the 
requirements for this sequencing is given in reference 2. 

The control task for the attitude control system of the test was to change 
the orientation of the model from the initially vertical position to a predeter- 
mined orientation defined by the reference angles 9g, \|rg, and The orien- 

tation to be achieved was established by precessing the gyro spin axes so that 
after precession these spin axes defined an orthogonal triad, angularly dis- 
placed from the x E~-» Y E-> and Z-g-axis system by the angles 0g, \|rg, and 

(See fig. 6.) The gyro precession event occurred at t = seconds. (See 

table.) The numerical values for the reference angles for the test were 
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®R = -9*5° (pitch down) 
t R = ^.7° (yaw right) 

0 R = -2.5° (roll counterclockwise) 

Upon achievement of the desired inertial reference angles, the responsi- 
bility of the attitude control system was to maintain this orientation subject 
to disturbances. 

The disturbances consisted of intentional moments acting independently 
about each rotational axis. Disturbance rockets were located on the dynamic 
model so that one produced a pitch-up disturbance, one produced a yaw- left dis- 
turbance, and the two remaining rockets produced a couple that gave a clockwise 
roll disturbance. The location as well as the thrust level of the rockets was 
such that pitch and yaw moments of about 2^7 foot-pounds and a roll moment of 
about 26 foot-pounds were realized. The magnitude of the disturbance moments 
was consistent with the 3 (7 critical disturbance characteristics of the space- 
vehicle upper-stage rocket motor during flight operation. The disturbances 
were introduced on a planar basis with the yaw disturbance first, followed by 
the roll disturbance, and then the pitch disturbance. 


ANALYSIS AND RESULTS 


General 

The results of the dynamic-test- stand investigation were generated from 
photographic data and telemetry data. The photographic-data reduction proce- 
dure is somewhat involved and a complete description follows. The telemetry- 
data reduction was accomplished in a more conventional manner. The results 
obtained from each data source will be given in this section. 


Camera Results 

A camera was mounted on a structure erected around the test stand so that 
it viewed the model from above. Placed between the camera and the model was a 
fixed dimensionally graduated framework. This framework represented the 
y E' Z E plane of the x E -> y e-> and Zg-axis system. The origin of the earth- 

fixed axis system lies on a fixed vertical which passes through the point of 
rotation of the dynamic model (bearing point). (See fig. 12.) 

The nose section of the spacecraft had two black lines painted on it which 
crossed at a point on the nose that corresponded to a point along the longitu- 
dinal axis of the dynamic model. The two painted lines were in mutually per- 
pendicular planes, and a black ball was fastened to the nose of the spacecraft 
at the point of intersection of the two lines. Figure 7 shows a typical photo- 
graph taken by the camera. 
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The purpose of the camera was to monitor the motion of the position- 
indicating hall so that data could he measured which would define the roll, 
pitch, and yaw Euler angles relative to the earth-fixed axis system. (See 
fig. 6.) The photographic data needed to define the pitch and yaw Euler angles 
are the coordinates of the hall in the y E> z E plane. The rotational orienta- 
tion of the hlack stripes in a plane perpendicular to the line of sight of the 
camera is necessary to determine the roll Euler angle. The equations required 
to determine the Euler angles 8, i|r, and (f) from the photographic data are 
derived and presented in appendix A as equations (Al8), (A19), and (A20) . 

The photographic data were also employed in the determination of attitude 
error angles. (See fig. 6.) These angles define the orientation of the dynamic 
model relative to the reference axis system. The reference axis system is known 
relative to the earth-fixed axis system hy the angles 6p, \)rp, and and the 

orientation of the dynamic model is defined relative to the earth-fixed axis 
system hy the photographically determined Euler angles. Therefore, the attitude 
error angles can be determined with use of the calculated Euler angles and 
defined reference angles . Equations for the attitude error angles 0 e , i|r e , 

and f & are derived and presented as equations (B7) , (B8) , and (B9) in 
appendix B. 

The results that were obtained from the photographic data are given in 
figure 8. These results were obtained by processing data at 0.125-second 
intervals, which corresponds to using every eighth frame obtained by the cam- 
era. Roll results are presented only during roll control because of inaccu- 
racies in the roll data after three-axis control when the nose of the model is 
rotated away from the line of sight of the camera. Figure 8(a) shows the Euler 
angles that were determined by using equations (Al8), (A19), and (A20). This 
plot gives an attitude time history of the dynamic model during the experi- 
mental test. The zero value for each parameter refers to the initial orienta- 
tion, and the dashed lines represent the inertial reference angles which define 
the attitude to be achieved. 

The closed- loop system comprised of the research control system, guidance 
system, and model exhibits type-0 servosystem steady- state characteristics 
(ref. 3)j that is, in order to produce a steady control deflection to balance 
a steady external moment, a steady actuating (attitude) error must be present. 
For this reason an attitude error is present during the time period the dis- 
turbance rockets are functioning. This error is represented by the deviations 
of the Euler angles from the reference values and occurs in each plane corre- 
sponding to the presence of a planned disturbance. (See fig. 8(a).) 

Figure 8(b) shows the variation of the roll, pitch, and yaw attitude error 
angles with time. These angles were computed from camera data with the use of 
equations (B7)., (B8), and (B9). Zero values for the gyro angles correspond to 
the dynamic model alined with the reference axis system. The steady- state 
error discussed previously is shown in figure 8(b) . 
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Telemetry Results 


Control-rocket thrust .- Pressure transducers were utilized to monitor the 
chamber pressure of each control rocket , and the results were telemetered to a 
receiving station. The chamber pressure was utilized to compute the thrust or 
available control force of the control rockets. The thrust was determined by 
using equations given in reference 4. The transducer of control rocket number 2 
did not function properly; therefore, only three thrust time histories were 
determined and these are shown in figure 9» The dashed curve represents an 
average of the experimentally determined results. 

Control parameters .- Figure 10 presents time histories of the attitude 
error angles, body rotational velocity components, and control-rocket positions 
that were obtained by telemetry. The attitude error angles are given in fig- 
ure 10(a). The outputs of the pitch and yaw attitude gyros were grounded to 
zero during roll control, which was necessary since the pitch and yaw attitude 
information is not used for control during this period. 

Figure 10(b) gives the variation with time of the body rotational velocity 
components p, q, and r which were sensed by rate gyros contained in the 
guidance system. The transients associated with release, roll capture, and 
three-axis capture can be seen in this figure. 

Figure 10(c) shows the position variation of each control rocket with time. 
These data were obtained by telemetry from the position-indicating potentiom- 
eters geared to the control rockets. For each control rocket one deck of the 
potentiometer provided data within ±10° boundaries to give good position-data 
resolution. The other data deck measured positions between ±60°. These two 
sources were combined to give the results presented in figure 10(c) where, 
for control rocket number 4, no data are shown beyond the ±10° boundaries. The 
steady-state deflections needed to control the intentional disturbances can 
also be seen in figure 10(c). 


COMPARISON AND DISCUSSION OF ANALOG AND EXPERIMENTAL RESULTS 


A graphical comparison between the results obtained from the experimental 
test and the results obtained from the analog simulation is given in figure 11 
where the experimental results are a repeat of previous presentations in this 
report. The analog simulation makes use of the equations presented in 
appendix C. 

Figure ll(a) shows a comparison of the Euler angles determined photograph- 
ically (experimental) with those determined by the analog analysis. An exami- 
nation of the photographically determined roll Euler angle shows that the model 
rolled positively during the time interval between release and the beginning of 
control. This positive roll was apparently the result of an extraneous torque, 
possibly produced by the separation mechanism. This motion is not shown in the 
analog results because the moment was not defined and was not included in the 
analog simulation. 
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Time histories of the pitch Euler angle show that the model was under the 
influence of a negative pitching moment, which was probably the result of winds 
and other external disturbances. This extraneous moment was estimated to be 
about -55 foot-pounds and was included in the analog simulation. The camera- 
determined yaw Euler angle diverges (t « 1 PP sec) from the defined reference 
value represented by the dashed line. This divergence is apparently caused by 
a drift of the gyro spin axis which defines the inertial reference, (in order 
to represent this drift in figure 8(a), the dashed line would be curved upward.) 
Since a constant reference was used in the analog analysis, this divergence is 
not noted in the analog results. 

The analog and camera Euler-angle results agree extremely well except for 
the explainable differences which are caused primarily by detailed disturbances 
that are not definable in the analog representation. The primary criteria for 
comparison are dynamic stability and damping characteristics. 

Figures ll(b), ll(c), and ll(d) give the variations of the attitude error 
angles with time. Since the experimental results were obtained from both telem- 
etry and photography, a three-way comparison was possible between these results 
and those from the analog analysis. Because the telemetry results are unavail- 
able during roll control, a portion of the comparison is lost. Much of the 
same basic discussion presented with the Euler angle results applies to the 
error-angle comparison, and as with the Euler angle comparison, the agreement 
is good. Although there are numerical differences, especially with the telem- 
etry data, the damping characteristics obtained from all three sources are 
comparable. 

Figure ll(e) gives a comparison of the rotational velocity components p, 
q, and r obtained by telemetry (experimental) and those obtained by analog 
simulation. Two pronounced points of disagreement are at roll capture and pitch 
capture. The roll capture disagreement was previously explained as a disturb- 
ance which was not included in the analog simulation. The same explanation 
applies to the pitch-capture disagreement as exemplified by the fact that the 
rate q was returning to zero when a disturbance redirected it in a negative 
direction. Much of the low-amplitude oscillation noted in the figure results 
from the control rockets deflecting indeterminately within unintentional and 
uncorrectable deadbands. A detailed comparison of the rotational -velocity 
results is not considered to be realistic; however, a survey of all transients 
shows good dynamic -stability agreement. 

A comparison of the control-rocket position time histories is given in 
figures ll(f) and ll(g). The oscillatory wave which appears to be riding on 
the deflection traces is the result of system nonlinearities. The steady-state 
deflections needed to control the intentional disturbances are seen by reviewing 
the figures. The graphical comparison of the control- rocket deflection results 
shows good agreement. 
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CONCLUDING REMARKS 


The results obtained from a dynamic-test- stand investigation of a space- 
vehicle upper stage employing actual electro-mechanical systems have been pre- 
sented. Experimental data received from an onboard telemetry system and a 
unique photographic technique were presented and compared with the results 
obtained from an analog s imul ation of the dynamical system. 

The Euler angles determined from photographic data showed good agreement 
with the analog- study results as did the attitude error angles computed from 
the camera data and those obtained in the analog study. A comparison between 
the attitude error angles obtained from these data sources and those reduced 
from telemetry data also showed good agreement in most areas. The body rota- 
tional velocity components obtained from telemetry and those obtained from the 
analog simulation showed the same general trends. The control- rocket position 
histories from these two data sources were in reasonable agreement. These com- 
parisons also showed that although telemetry results were somewhat limited in 
accuracy, the telemeter data can be used for a comprehensive analysis of the 
dynamical system. 

The analog simulation employed herein was also used to analyze and define 
system gains and sequences in a previous analysis of the space-vehicle upper 
stage. Since this dynamical system functioned as predicted by the analog study, 
it is concluded that the analog simulation is an adequate mathematical repre- 
sentation of the dynamical system. 


Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., May 12, 19^5 • 
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APPENDIX A 


EQUATIONS FOR REDUCING CAMERA DATA 


Equations are developed in this appendix which define the attitude of the 
dynamic model as a function of parameters obtained from photographic data. The 
data obtained by the camera technique enployed in the dynamic-test-stand inves- 
tigation are the coordinates of P* (see fig. 12) and a protractor-determined 
roll angle. The coordinates of P' do not represent the true coordinates of 
the position- indicating ball in the Yg,Zrg plane due to parallax. The true 

coordinates are represented by the point P. The measured roll angle also does 
not represent the true roll angle , because the line of sight of the camera is 
not along the longitudinal axis of the dynamic model. Because the raw photo- 
graphic data must be corrected prior to determination of the desired attitude 
angles, the correction equations will be presented first followed by a presen- 
tation of the attitude equations. 


DATA CORRECTION EQUATIONS 


The Initial problem is one of defining relationships so that the coor- 
dinates of P 1 can be corrected to give the coordinates of P, and the 
protractor-measured roll angle can be corrected to give the true roll angle. 

The two right triangles shown In figure 12 with the Yg-axis as one leg, 
apexes at point 0, and hypotenuses OP’ and OP are similar j therefore. 

Ye _ Ye 

OP 7 OP 

and 

Z E _ fE_ 

OP 7 OP 

These equations can be modified to give the following: 

■ OP (A3) 

to) 2 



(Al) 


(A2) 
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and 

Zrp = OP (A4) 

+ ( z e ) 2 

In order to solve equations (A3) and (A^) It is necessary to define OP. 
From figure 12 it can be seen that 

OP = r m sin p (A5) 

where r m is known and a relationship for sin p must be determined. 

Equation (A 6) is the law of sines for the triangle containing the angles 
a, p, and co. 


r m = CB = h + e + r m ( A 6) 

sin a sin P sin co 


From equation (A 6) 


sin p 


CB 


sin a 


(AT) 


and CB can be defined as 

CB = (h + e + rin)cos a + r m cos co (A8) 

Also, from equation (A6) 


h + e + r m 

sin co = — sin a 

x m 


and from figure 12 sin a may be expressed as follows: 

OP 7 


(A9) 


sm a = 


7 


(A10) 


OP 1 


Use of equations (A9) and (A10) in substituting for the cosine functions of 
a and co in equation (a 8) gives 


CB 


- ( 


h + e + r. 


my 


f _ 

i _ p\ 


9 P 


(h 2 + OP’ 2 ; 

1 - ( h + e + r m ) 

rop* 


(All) 
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The negative sign of the second radical of equation (All) is necessary since 
cd is greater than 90° for all applications in this report. Combining equa- 
tions (A7), (AlO), and (All) gives 


sin p = 


0P T 


r m (h 2 + OF 2 ) 


h(h + e + r m ) 


- jr m 2 (h 2 + OP’ 2 ) - 0P' 2 (h + e + r m ) 2 


(A12) 


Equations (A3), (A4-), (A5), and (A12) can now he combined to give 


y E = y E\ 


h(h + e + r m ) - | 

r 2 
• L m 

h2 + ( y E)' 

' + 

[ z e) 2 

- 

* (•if 

(h + e + r m f 

V 

h 2 + 


) 2 + 

K ) 2 

J 


(A13) 


and 


h(h + e + r m ) - )Jr m 2 h 2 + (y E )‘ 

3 

' + 

[■if 

- 

( y Ef + ( z e) 2 ( h + e + r m f 

h 2 + 

( y il 

f + 

( z e) 

J 

(Al4) 


Equations (A13) and (Al4) define the true coordinates y„ and z„ of the 

ill 

position-indicating ball in the Y E ,Z E plane. 

The protractor-measured roll angle 0* is in a plane perpendicular to the 
line of sight of the camera. The true roll angle 0 is in a plane perpendic- 
ular to the longitudinal axis of the dynamic model. Therefore, the measured 
roll angle must be projected into the desired plane. The following equation 
defines this projection: 

sin 0* = cos T} sin f (A15) 

The value of cos T] is the unknown in equation (A15) which must be defined. 
Since T) and cd are supplementary angles, tj may be expressed in terms of 
cu as 

T| = 180 - CD (Al6) 


Ik 
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and by combining equations (Al6) and (A9) 


cos T) = 


E? 

ru 

(h 2 + OP' 2 ) 

- ( h + e 

+ r m ) 2 0P' 2 

1 r 2 I 

± m 

(h 2 + 0 F 2 ) 

1 


(A17) 


Equations (A15) and (A17) can be used to define the true roll angle 0 as 
follows : 


sin 0 = sin 0' 


1 ^ 

j ■ L m 

h2 + (ye 

f + frf. 


r 2 
j x m 

» 2 - {vkf * (■ 

1 \ 2 ~ 

z e) 

-( 

h + e + r m) 2 

( y E) + ( z e) 


(A18) 


Equations (A13), (A14), and (Al8) yield the results needed to define the atti- 
tude of the dynamic model as a function of the raw photographic data. 


ANGULAR AOTITUDE EQUATIONS 


Euler angles 0, \|r, and 0 define the attitude of the dynamic model in 

the X E -, Y-g-, and Z E -axis system. The sine of the roll Euler angle was given 

by equation (Al8), and the pitch and yaw Euler angles are defined by the fol- 
lowing equations: 


0 



(A19) 


* = 



(A20) 
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APPENDIX B 


EQUATIONS USED IN DETERMINING ATTITUDE ERROR ANGLE 


Equations were derived in reference 2 which defined the attitude error 
angles in terms of Euler angles and a pitch reference angle. A similar deriva- 
tion is presented in this appendix in order to define the attitude error angles 
in terms of Euler angles and roll, pitch, and yaw reference angles. The deriva- 
tion consists of a series of orthogonal transformations between axis systems. 
Figure 6 serves as a useful tool in reviewing the following derivation. 

The transformation between the earth-fixed axis system ^Xg, Yg, and Zg^ 
and the body axis system is 


X 


1 

0 

0 

— 

cos 0 

0 

-sin 0 


COS ijf 

sin T|r 

0 


X E 

Y 

= 

0 

cos 0 

sin 0 

0 

1 

0 


-sin i|r 

cos i|r 

0 


y e 

Z 


0 

-sin 0 

cos 0 

sin 0 

0 

cos 0 


0 

0 

1 


Zg 


(Bl) 


which, when simplified, becomes 


X 


cos 0 cos \|r 

Y 


sin 0 sin 0 cos i|r 



- cos 0 sin i|r 

2 


cos 0 sin 0 cos \jr 



+ sin 0 sin \|r 


cos 0 sin i|f 

sin 0 sin 0 sin t|t 
+ cos 0 cos f 

cos 0 sin 0 sin ^ 
- sin 0 cos \| r 


-sin 0 


V 

sin 0 cos 0 


y e 

cos 0 cos 0 


Z E 

_ 


_ — 


(B2) 


The transformation between the earth-fixed axis system and the reference 
axis system is 


X R 


COS 0g COS l|fg 

cos 0g sin ij/g 

Yr 

= 

sin 0g sin 0g cos i|r R 
- cos 0p sin i|rg 

sin 0g sin 0g sin \|fg 
+ cos 0p COS 

Z R 


cos 0g sin 0g cos \|rg 
+ sin 0g sin \|rg 

cos 0g sin 8g sin i|rg 
- sin 0g cos i|rg 


-sin 9g 


X E 

sin 0g cos 0g 


y e 

COS 0g COS 0g 


Z E 

j 




(B5) 
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By transposing equation (Bjj) and substituting the 
tion (B2) , the following equation is obtained: 


results into equa- 


X 


z 


COS 8 

cos if 

COS 0 R COS 


COS 0 

cos \|r 

sin 0 R sin 0 R cos i^ R 


cos 0 

COS i|r 

cos 0 R sin 0 R cos \jr R 



+ COS 

0 sin 

cos 

e R sin *R 


- 

COS 

9 cos 

* 

cos 

0R aln +R 


+ 

cos 

0 cos 


sin 0 R sin >Jr R 



+ sin 

0 sin 

e R 



+ 

cos 

0 sin 

* 

sin 

0r sin 0 R sin i^r 


+ 

cos 

0 sin 

* 

cos 0 R sin 0 R sin \jf R 


X R 






+ 

cos 

9 sin 

if 

cos 

0R COS l|f R 


- 

cos 

0 sin 

if 

sin 0 R cos iJr R 








- 

sin 

0 sin 

h 

cos 0 R 



sin 

0 cos 

0R 0 R 



sin 0 

sin 0 

COS if 

cos 0 R cos 

♦r 

sin 0 

sin 0 

cos if 

Bin 0 R sin 0 R cos 

*R 

sin 0 

sin 0 

cos cos 0 R sin 0 R cos 

*R 


- cos 

0 sin 

cos 

e R cos + R 


- 

sin 

0 sin 

0 

cos 

y cos 0 R sin \Jr R 


+ 

sin 

0 sin 

0 

cos if sin 0 R sin i|r R 



+ sin 

0 sin 

0 sin 

i | t cos 0 R sin if R 

- 

cos 

0 sin 

if 

sin 

0 R sin 0 R cos \|f R 


- 

cos 

0 sin 

if 

cos 0 R sin 0 R cos iJr R 



+ cos 

0 cos 

if cos 

e R sin i|r R 


+ 

cos 

0 sin 

if 

cos 

0R = ln tR 


- 

cos 

0 sin 

if 

sin 0 R sin ty R 



- sin 

0 cos 

G sin 

0 R 


+ 

sin 

0 sin 

0 

sin 

if sin 0 R sin 0 R sin ty R 

+ 

sin 

0 sin 

0 

sin if cos 0 R sin 0 R sin \|r R 

y r 






+ 

sin 

0 sin 

0 

sin 

cos 0jj cos \|r R 


- 

sin 

0 sin 

0 

sin i]r sin 0 R cos ijr R 








+ 

cos 

0 cos 

if 

sin 

0r sin 0 R sin \jr R 


+ 

cos 

0 cos 

if 

cos 0 R sin 0 R sin >Vr 








+ 

cos 

0 cos 

if 

cos 

0R <=° B 


- 

cos 

0 cos 

* 

sin 0 R cos \|r R 








+ 

sin 

0 cos 

0 

sin 

0R COB e R 


+ 

sin 

0 cos 

0 

cos 0 R cos 0 R 



cos 0 

sin 0 

cos 

cos 0 R cos 

+ R 

cos 0 

sin 0 

COB \|r 

sin 0 R sin 0 R cos 

♦r 

cos 0 

sin 0 

cos if cos 0 R sin 0 R cos 

♦s 


+ sin 

0 sin 

\|r cos 

0 R cos \|r B 


- 

cos 

0 sin 

0 

cos 

cos 0 R sin i|r R 


+ 

C08 

0 sin 

0 

cos \|r sin 0 R sin ty R 



+ cos 

0 sin 

0 sin 

t cos 0 R sin if R 

+ 

sin 

0 sin 

if 

sin 

0r sin 0 R cos 


+ 

sin 

0 sin 

if 

cos 0 R sin 0 R cos \{r R 



- sin 

0 cos 

cos 

0 R sin \|/ R 


- 

sin 

0 sin 

if 

cos 

0R s ln +r 


+ 

sin 

0 sin 

if 

sin 0 R sin \|r R 



- cos 

0 COB 

0 sin 

e R 


+ 

cos 

0 sin 

0 

sin 

>|r sin 0 R sin 0 R sin \|f R 

+ 

cos 

0 sin 

0 

sin cos 0 R sin 0 R sin \|r R 

Z R 






+ 

cos 

0 sin 

0 

sin 

>Jf cos 0j^ cos ^r R 


- 

cos 

0 sin 

0 

sin \|r sin 0 R cos \|r R 








- 

sin 

0 cos 

if 

sin 

0 R sin 0 R sin 


- 

sin 

0 cos 

if 

cos 0 R sin 0 R sin y R 
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sin 

0 cos 

* 

cos 

0R =° s *R 


+ 

sin 

0 cos 

V 

sin 0 R cos \Jr R 








+ 

cos 

0 cos 

0 

sin 

0R e R 


+ 

cos 

0 cos 

0 

cos 0 R cos 0 R 

_ 



(Bb) 


Equation (B^) represents a transformation between the body axes and the refer- 
ence axes. An additional transformation between the body axes and the refer- 
ence axes is as follows: 


X 


cos 0 


£ cos i|r* 


-sin 9, 


sin 

0 € sin 

0 £ cos 

ife 

sin 

0 £ sin 

0 £ sin 

ift 

sin 0 e cos 0 £ 

- 

cos 0 e 

sm \(r £ 


+ 

cos 0 £ 

cos \|/ £ 



cos 

0 e sin 

0 € cos 

** 

cos 

0 e sin 

0 e sin 

i * 

cos 0 £ cos 0 £ 

+ 

sin 0 e 

sin i|r* 


- 

sin 0 £ 

cos 




% 


R 


J R 


(B5) 


Since the attitude error angle is measured in the plane of the dynamic 
model and the yaw angle i|r £ noted in the transformation is the projection of 

the attitude error angle in the reference X rAr plane, the following expres- 
sion is necessary to define the yaw attitude error angle: 

sin \|r £ = sin \j/ £ cos 0 £ (B6) 
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By equating the equal elements of (B^) and (B5), the following attitude-error- 
angle equations are obtained: 

sin 0 e = -cos 0 cos ijr cos 0 R sin 0 R cos i|r R 

- cos 0 cos \|r sin 0 R sin i|r R 

- cos 0 sin •vjr cos 0 R sin 0 R sin i|r R 

+ cos 0 sin i|r sin 0 R cos i|r R 

+ sin 0 cos 0 R cos 0 R (D7) 

sin = cos 0 cos i|r sin 0 R sin 0 R cos i|/ R 

- cos 0 cos i|r cos 0 R sin i|r R 

+ cos 0 sin \|r sin 0 R sin 0 R sin * R 

+ cos 0 sin Ajr cos 0 R cos i|r R 

- sin 0 sin 0 R cos 0 R (®8) 

sin 0 6 cos 0 € = sin 0 sin 0 cos i|r cos 0 R sin 0 R cos t R 

+ sin 0 sin 0 cos i | r sin 0 R sin i|r R 

- cos 0 sin T|r cos 0 R sin 0 R cos l|f R 

- cos 0 sin ijr sin 0 R sin i]r R 

+ sin 0 sin 0 sin i|r cos 0 R sin 0 R sin i|r R 

- sin 0 sin 0 sin i|r sin 0 R cos i|r R 

+ cos 0 cos i]r cos 0 R sin 0 R sin i)r R 

- cos 0 cos i|r sin 0 R cos i|r R 

+ sin 0 cos 0 cos 0 R cos 0 R (B9) 
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APPENDIX C 


EQUATIONS USED IN ANALOG COMPUTER ANALYSIS 


The mathematical descriptions of the dynamical systems used to generate 
the results presented in this report are the same as those given in reference 2 
except for the vehicle and trajectory simulation. The representation for the 
vehicle and trajectory simulation of this investigation is given by the fol- 
lowing equations and differs from that of reference 2 because no translation 
occurred during the dynamic test: 


Six 

= -F ^sin 

- sin - sin &2 + sin 5]^Xp + My 

(Cl) 

qiy + 

rplx - rpl z 

= -F^sin &2 + s l n ^4)^ x cr " x q) + My 

(C2) 

rl z + 

pqI Y - pql x 

= -F^sin + sin Sj^x cr - x r ^ + Mg 

(C3) 



q sin 0 + r cos 0 
cos 0 

(c4) 



0 = p + ijr sin 0 

(C5) 


• 

0 

= q cos 0 - r sin 0 

(c 6) 



t = / t dt + t Q 

J t 

(C7) 



0 = / 0 dt + 0 O 

d t 

(C8) 



0 = J 0 dt + 0 O 

(C9) 


19 


REFERENCES 


1. Fall, R. C.; Koch, Robert L.j and DeBord, C. J.: Development of the Rocket 

Velocity and Attitude Control System. Eng. Dept. Rept. No. 2137 (Contract 
No. NAS 5-433) , Allison Div., Gen. Motors Corp., May 29, 196l. 

2. Young, A. Thomas; and Harris, Jack E.: An Analog Study of a Rotating-Solid- 

Rocket Control System and Its Application to Attitude Control, of a Space- 
Vehicle Upper Stage. NASA TN D-2366, 1964. 

3. Nixon, Floyd E.: Principles of Automatic Controls. Prentice-Hall, Inc., 

C.1953- 

4. Barrere, Marcel; Jaumotte, Andre; De Veubeke, Baudouin Fraeijs; and 

Vandenkerckhove , Jean: Rocket Propulsion. Elsevier Pub. Co., i960. 


1 


1 


20 










Figure 3.- Block diagram of major elements of dynamic model. 








Top view 

(looking in +Z direction) 


Side view 


(looking i 


+Y direction) 


Figure 5.- Sketch showing control rocket deflection convention. 





Figure 7.- Typical view from top camera. Each graduation equal to 3 inches. 
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(a) Variation of 0, e, and i/i with time. 

Figure 11.- Comparison of analog simulation results and dynamic-test-stand investigation results. 
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(b) Variation of fi £ with time. 



(c) Variation of with time. 
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Figure 11.- Continued. 



(e) Variation of p, q, and r with time. 
Figure 11.- Continued. 
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